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SUMMARY 


A classification of all the compensators that stabilize a given general plant in a linear, time-invariant 
multi-input, multi-output feedback system is developed. This classification, along with the associated 
necessary and sufficient conditions for stability of the feedback system, is achieved through the intro- 
duction of a new parametrization, referred to as R-Parametrization, which is dual of the familiar 
Q-Parametrization. The classification is made according to the stability conditions of the compensators and 
the plant by themselves; and the necessary and sufficient conditions are based on the stability of Q and R 
themselves. 


1. INTRODUCTION 


Recently, there is a new formulation of the old problem of designing a compensator for a given 
time-invariant multi-input, multi-output (MIMO) plant to meet a set of performance specifications. This 
new formulation (e.g., Callier et al. (1982), Midyasagar (1985), and Boyd et al. (1988)) takes the view of 
choosing the compensator from among all die compensators that stabilize the given plant and therefore 
represents a significant departure from the earlier design methods such as Linear-Quadratic-Gaussian 
(LQG) (Bryson et al. (1969)), Inverse Nyquist Array (Rosenbrock (1969)), and Characteristic Loci 
(MacFarlane et al. (1977)). Essential to this new formulation is a deeper understanding of the nature and 
the characteristics of all die compensators that stabilize any given plant. This paper is therefore addressing 
this specific area. 

Consider a multivariable unity feedback system as shown in figure 1. The dimension of the com- 
mand input ul (and the plant output y2), and the disturbance input u2 (and the compensator output y 1) 
are nO and ni respectively. It is assumed that (i) the compensator, C(s), and the plant, P(s), are in gen- 
eral proper, (ii) neither the plant nor the compensator has unstable hidden modes, and (iii) the system is 
well-posed (i.e., det(I + DpDc) * 0). The transfer matrices from utoy and from utoe are given 
respectively in equations (1) and (2). 


*This work was performed during the academic year, 1987-88 while the author was with the Aeronautics/Astronautics 
Department., Stanford University, as a visiting professor. 


1 



COMPENSATOR 
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u2 



Figure 1- Multivariable unity feedback system. 


C(I + PC)- 1 

-CP(I + CP)- 1 

PC(I + PC)- 1 

P(I + CP)- 1 

(I + PC)" 1 

-P(I + CP)- 1 

_ C(I + PC)- 1 

(I + CP)- 1 


( 1 ) 

( 2 ) 


It is known (e.g., Callier et al. (1982)) that these transfer matrices can be cast into a form which is affine 
in Q as shown in equations (3) and (4). This "linearizing" process is accomplished using the nonlinear 
transformation pair, equations (5) and (6), and is commonly referred to as Q-Parametrization. 1 


H yu = 

Q 

-QP 

PQ 

P(I - QP) J 

Heu = 

' I-PQ 
Q 

-P(I-QP)' 

I-QP 


(3) 

(4) 


Q = C(I + PC)" 1 = (I + CP)” 1 C 


(5) 


C = Q(I - PQ)" 1 = (I - QP) _1 Q (6) 

It is worth noting, from equations (l)-(4), that Q is the transfer matrix from ul to yl (or from ul to e2), 
i.e., 


Q — Hyl )U l - H^ul (2) 

We will defer the discussion of other significant characteristics of the Q-Parametrization to a later section 
in this paper. For now we will proceed to introduce another parametrization, a new one which will hence- 
forth be referred to as R-Parametrization. 


^The idea of the Q-Parametrization for SISO systems was discussed in the literature many years ago (e.g., Newton et al. 
( 1957 )). 
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2. R-PARAMETRIZATION 


This linearizing process was motivated by the problem of answering the question: For a given 
stable C(s) (i.e., C(s) by itself is stable), what are all the plants, P, stable or otherwise, which are 
stabilizable by C(s)? This is "dual" to the problem of answering the question: Given a stable P, what are 
all the compensators, C, that can stabilize P? (Note that C's may or may not be stable by themselves.) 

We note that this latter question has been answered by the Q-Parametrization (Zames (1981), Desoer et al. 
(1981)) described earlier. 

To seek a linearizing process directly aiming at answering the first question, it is proposed that the 
following nonlinear transformation pair, equations (8) and (9), be used: 

R = P(I + CP)" 1 = (I + PC)" 1 P (8) 

P = R(I - CR) _1 = (I - RC) -1 R (9) 

It is important to recognize that R is the transfer matrix relating the plant output to disturbance input u2. 
With this pair of transformations, it can be shown that the transfer matrices (1) and (2) become: 


H yu = 
Heu = 


C(I-RC) 

RC 

I-RC 

C(I-RC) 


-CR 

R 

-R 

I-CR 


( 10 ) 

(ID 


It is seen from equations (10) and (11) that the closed-loop transfer matrices are now affine in R. From 
equations (3), (4), (10), and (11), the following identities are evident: 


R = P(I - QP) = (I - PQ)P 

(12) 

Q = C(I - RC) = (I - CR)C 

(13) 

PQ = RC 

(14) 

QP = CR 

(15) 


These identities will prove useful in establishing some of the important stability criteria for the feedback 
system to be discussed next. 


3. STABILITY CRITERIA FOR THE FEEDBACK SYSTEM 


Since we are concerned with compensators and plants which are proper (not just strictly proper), 
thereby broadening the engineering applicability of our results, we will be concerned in this paper with 
bounded-input-bounded-output (BIBO) stability. The following fact is known. 

FACT1: If the plant P is stable by itself, then C stabilizes P (i.e., the feedback system is stable in the 
sense of BIBO) iff Q is stable (e.g., Callier et al. (1982)). 
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Fact 1 is clear from equation (3) or (4). An engineering interpretation for that seemingly mysterious nec- 
essary and sufficient condition can be made by simply redrawing figure 1 in the form of figure 2. Since 
the given plant P is stable, d is bounded for bounded u2. We note therefore that the signal d is equiv- 
alent to -ul. Thus, all the signals in the system are bounded iff Q is stable (recall from equation (7) that 
Q is the closed-loop transfer matrix from command input ul to the compensator output yl). The fol- 
lowing fact is "dual" of Fact 1. 



Figure 2 - Multivariable unity feedback system in terms of P and Q. 

FACT 2: If the compensator C is stable, then C stabilizes P iff R is stable. 

This fact is clear from equations (10) and (11). A somewhat less visualizable proof of this fact is to show 
algebraically from equation (1) or (2) that the feedback system is stable iff P(I + CP) -1 is stable as is 
usually done in the literature (e.g., Vidyasagar (1985)). Again, a good physical insight into this fact can 
be made by redrawing figure 1 in the form of figure 3. 



Figure 3.- Multivariable unity feedback system in terms of C and R. 

We note that Cul is a bounded input similar to u2 since C is stable, and that the closed- loop transfer 
matrix from u2 to y2 is simply equal to R. 

FACT 3: If both C and P are stable by themselves, then C stabilizes P iff (I + PC) -1 (or alternatively, 
(I + CP) -1 ) is stable. 
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The necessity is clear from equation (2). For sufficiency, we see that (1) can be arranged in the forms of 
equations (16) and (17) and the fact is established. 



C(I + PC)" 1 
PC(I + PC)” 1 


-C(I + PC)" 1 P 
(I + PC)" 1 ? 


(16) 



(I + CP )" 1 C 
P(I + CP )- 1 c 


-CP(I + CP)- 1 
P(I + CP)- 1 


(17) 


From equations (2), (4), and (1 1) we see that (I + PC) -1 = I — PQ = I - RC. This, along with equa- 
tions (5) and (8), indicates that (I + PC) -1 is stable iff either Q or R is stable, since both C and P are 
stable. We have, therefore, 

FACT 3A: If C and P are both stable, then C stabilizes P iff either Q or R is stable. 


We next establish Fact 4. 

FACT 4: C and P are both unstable, but they do not have common unstable poles. C stabilizes P iff 
both QandR are stable. 


Necessity is clear by combining equations (4) and (1 1) in the following form 


Heu = 


I-PQ 

-R 


I-RC 

-R 

. Q 

I-QP . 


. Q 

I-CR . 


or in the form of equation (19), 


Q 

-QP 


Q 

-CR 

PQ 

R 


RC 

R 


(18) 


(19) 


For sufficiency, we see from the first equality of equation (18) or (19) that if instability occurs in the feed- 
back system it must be caused by the unstable poles of P, since both Q and R are stable. However, the 
right-hand side of equation (18) or (19) indicates that the instability can only be from the unstable poles of 
C. Since PandC have no common unstable poles, we conclude that the feedback system cannot be 
unstable if both Q and R are stable, thus establishing Fact 4. 


Fact 4 was first found by Desoer et al. (1975) using complex variable analysis, not in the context 
of Q- and R-Parametrization. To gain a better insight into this fact, we shall express the terms 
I-PQandl-QP in equation (18) in terms of only QandR. To do this, let X = I-PQ and 
Y = I - QP. Then from equations (9), (14), and (15), we have 

X = I - X -1 RQ (20) 


Y = I - QRY -1 


( 21 ) 


Solving these two second-order matrix equations yields 


x-I 

X “2 


I±(I-4RQ) 1/2 


( 22 ) 
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Y=i[l±(I-4QR) 


(23) 


l/2l 


One of the two solutions in equations (22) and (23) provides the answer for I - PQ and I - QP, 
respectively, as will become clear in the sequel. However, both solutions of X and Y are stable if 

(I-4RQ) 1 / 2 and (I-4QR) 1/2 

are stable. Note that (I - 4RQ) and (I - 4QR) are stable if both Q and R are stable. For the SISO case, 
we note that (1 - 4RQ) 1 /2 = (1 - 4QR) 1 / 2 and they are stable if (1 - 4RQ) is stable. This fact does 
not hold for the MIMO case, however. Indeed, from equations (12) and (13), we have 

I - 4RQ = (I - 2PQ) 2 = (I - 2RC) 2 
and 

I _ 4QR = (I - 2QP) 2 = (I - 2CR) 2 

We see that (I - 4RQ)V2 = I - 2PQ = I - 2RC or (I - 4QR)i/ 2 = I - 2QP = I - 2CR can be 
unstable if P and C have common unstable poles. Examples 1 to 4 in table 1 show the four combina- 
tions of the stability characteristics of (I - 4RQ) 1 ^ 2 and (I - 4QR) 1 #, given that Q and R are both 
stable. Example 1 shows that C and P do not have a common unstable pole; consequently, (I - 4RQ) 1 / 2 
and (I - 4QR) 1 / 2 are both stable. In examples 2 to 4, C and P have common unstable poles, which give 
rise to three combinations of instability in (I - 4RQ) 1 / 2 and (I - 4QR) 1 ^ 2 . It is interesting to observe 
from these examples that when the instabilities occur in (I - 4RQ) 1 / 2 or (I - 4QR)l/2, there are trans- 
mission zeros coinciding with the unstable common poles of C and P in the associated transfer matrices 
I - PQ, PQ, or I - QP, QP, respectively (see table 3). This property of having transmission zeros 
coinciding with poles, which has no parallel in a SISO system, may exist in a MIMO system, and is the 
underlying reason why (I - 4RQ) 1 / 2 or (I - 4QR) 1 '' 2 becomes unstable while both Q and R are stable. 

It may be appropriate at this point to return to equations (22) and (23), and examine the two 
solutions of X and Y. Let XI = [I + (I - 4RQ)!/ 2 ]/2 and X2 = [I - (I - 4RQ)!/ 2 ]/2; 

Y1 = [I + (I - 4QR)!/ 2 ]/2 and Y2 = [I - (I - 4QR) 1 / 2 ]/2. Then since (I - 4RQ)V2 = 1 _ 2PQ 
and (I - 4QR) 1 / 2 = I - 2QP, it follows that XI = I - PQ, Y1 = I - QP, and X2 = PQ, 

Y2 = QP. This further reinforces the fact that when (I - 4RQ) 1 / 2 is stable, both Heiul and H y 2 u i are 
stable, and when (I - 4QR) 1 / 2 is stable, He2u2 and H y 2 u l are stable (see eqs. (18) and (19)). 

It is important to point out that, given Q and R being both stable, (I - 4RQ) 1 / 2 and (I - 4QR) 1 / 2 
can sometimes remain stable (thus the feedback system remains stable) even if there are common unstable 
poles in C and P. In other words, when Q and R are both stable, having common unstable poles in C 
and P does not imply that (I- 4RQ) 1 / 2 or (I- 4QR) 1 / 2 must be unstable. This fact is illustrated by 
example #5 as listed in tables 1 and 2. Summarizing the above discussion, we can state that, given Q and 
R being both stable, and C and P being both unstable, (a) C always stabilizes P for the SISO case, 
and (b) C almost always stabilizes P for the MIMO case with the exception being for some "patho- 
logical" cases in which C and P have common unstable poles that happen to cause (I - 4RQ)^ or 
(I - 4QR) 1 / 2 to become unstable. Fact 5 in the following addresses more precisely those pathological 
cases. 

FACT 5: C and P are unstable and they have common unstable poles. C stabilizes P if (i) Q and R are 
stable, and (ii) (I - 4RQ) 1 / 2 and (I - 4QR) 1 / 2 are stable. 
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TABLE 1.- FIVE EXAMPLES OF MIMO UNITY FEEDBACK SYSTEMS. 








TABLE 1.- CONTINUED. 


1 1 _ 

s + 1 s 

„ 1 


0 - 
s 



b = s 2 + s+ l; c = s 2 + 3s+1 


s + 1 s(s + l)(s + 2) 
c be 



s(s + 1) s(s + l) 2 
c be 


(I-4RQ) 1/2 = 


b 2(s + l) 2 
c be 

« s 2 + s - 1 


(I-4QR) 1/2 = 


2(s + 1) 2 (2s + 1) 
sbe 

s 2 + s - 1 
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TABLE 1.- CONTINUED. 


Example 




b = s 2 + s + 1 ; c = s 2 + 3s + 1 




b 

c 

(I-4RQ) 1/2 = 

0 


2(s + l) 2 (2s+ 1) 
sbc 

s 2 + s - 1 
b 


b 2(s + 1) 3 

c be 

<I-4QR)W= 2 

0 S - 

b 

NOTE: This is the example 9 
of Desoer et al. (1975). 




TABLE 1.- CONTINUED. 



Example 


d = s 2 + 2s + 5 


■_2 1 _ 

s + 1 s- 1 
C= . 2 


2(s + 1) 
d 


2(s + 1) 
d 


P = C 

Note: This is the example 7 
of Desoer et al. (1975). 


(I-4RQ) 1/2 = 


(s + 3)(s - 1) 8(s + 1 y 

d d 2 (s - 1) 

(s + 3)(s - 1) 
d 


(I-4QR) 1/2 = (I-4RQ) 1 / 2 




TABLE 1.- CONCLUDED. 



Parameter 


#5: 


C = 


S 4- 1 
S 

0 


0 

2(s + 1) 


P = 


s+ 1 


0 - 
s 


Note: In this example, 
CandP have common 
unstable pole at s = 0. 


Q = 


e = s 2 + 2s + 2 


1 

0 


2(s + 1) 
e 

2s(s + 1) 


; R = 


(s + l) 2 
0 


(I-4RQ) 1/2 = 


(I-4QR) 1/2 = 


s- 1 
s + 1 

0 

s- 1 

S + 1 
0 


e(s + 1) 

s^ 

e 


4s 

e 

s 2 - 2s - 2 


_2s 

e 

s 2 - 2s - 2 


TABLE 2.- STABILITY CHARACTERISTICS OF Q AND R PARAMETERS 
FOR THE FIVE EXAMPLES IN TABLE 1. 


Example 


Parameter 

1 

2* 

3* 

4* 

5* 

Q 

stable 

stable 

stable 

stable 

stable 

R 

stable 

stable 

stable 

stable 

stable 

(I-4RQ) 1 / 2 

stable 

stable 

unstable 

unstable 

stable 

(I-4QR) 1 / 2 

stable 

unstable 

stable 

unstable 

stable 


*C and P have common unstable poles. 
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TABLE 3.- POLES AND TRANSMISSION ZEROS OF VARIOUS TRANSFER MATRICES FOR 

THE FIVE EXAMPLES IN TABLE 1. 


Example 

Parameter 

Poles 

Transmission zeros* 

#1 

(I-4RQ)l/2 

{-1; -(1/2)0 ±V7j)} 

{1; 0/2)(3 ±VI7)) 


I-PQ 

(-1; -(1/2X1 ±V7j )} 

{0; 0; 0) 


PQ 

{-1; -(1/2)0 ±V7j)} 

{-1} 


(I-4QR) 1 / 2 

{-1; -0/2)(l ±V7j)) 

{1; (1/2X3 ±VT7)) 


I-QP 

{-1; -(1/2)0 ±V7j)} 

{0; 0; 0} 


QP 

{-1; -0/2)(l ±V7j» 

{-1} 

#2 

(I-4RQ)V2 

{-(l/2)(3 ±V5); -0/2)(l ±V3j)} 

{-<1/2)0 ±V3j); -0/2)(l ±V5)} 


I-PQ 

{-0/2)(3 ±V5); -(1/2)(1 ±V3j)} 

{0; -1; -1; -1} 


PQ 

{-(1/2X3 ±V5); -(1/2)0 ± V3j)} 

{0} 


(I-4QR) 1 / 2 

{0; — <l/2)(3 ± V5); -(1/2)(1 ± V3j)) 

{0; -0/2)(l ± V3j); -(1/2)(1 ± V5» 


I-QP 

{0; — (l/2)(3 ± V5); -(1/2)(1 ± V3j)} 

{0; 0; -1; -1; -1} 


QP 

{0; -(1/2) (3 ± V5); -(1/2)(1 ± V3j)) 

(0; 0} 

#3 

(I-4RQ)V2 

{0; -<l/2)(3 ± V5); -{1/2)0 ± V3j)} 

{0; -(1/2)0 ± V3j); -(1/2)(1 ± V5)} 


I-PQ 

{0; -0/2)(3 ± V5); -(1/2)(1 ± V3j)} 

{0; 0; — 1; — 1; — 1 ) 


PQ 

{0; -(1/2) (3 ± ^); — (1/2)(1 ± V3j)) 

{0; 0} 


(I-4QR) 1 / 2 

{-(l/2)(3 ±V5); -0/2)(l ± V3j)} 

{-(1/2)(1 ±V3j); -0/2)(l ±V5)} 


I-QP 

{-(1/2X3 ±V5); -0/2)0 ±V3j)} 

{0; —1; —1; —1 } 


QP 

{-(l/2)(3 ±V5); -0/2)(l ±V3j)} 

{0} 


*The transmission zeros are defined in terms of the Smith-McMillan form (see, e.g., Sinha (1984)). 
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TABLE 3.- CONCLUDED. 


Example 

Parameter 

Poles 

Transmission zeros 

#4 

(I-4RQ)l/2 

{ 1; — 1 ±2j;-l ±2jJ} 

{1; 1; 1; —3; —3} 


I-PQ 

{l;-l±2j;-l±2j}} 

{ l; —l; —l; — l; — l } 


PQ 

{l;-l±2j;-l±2j}} 

{1} 


(I-4QR)l/2 

{ 1; — 1 ±2j;-l ±2j}} 

{1; 1; 1; —3; —3} 


I-QP 

{ 1; — 1 ±2j;-l ±2j}} 

{ l; —l; —l; — l; — l } 


QP 

{ 1; — 1 ±2j;-l ±2j}} 

m 

#5 

(I-4RQ) 1 / 2 

(—1; —l ± j) 

U;1±V3} 


I-PQ 

{— i; — i ±j) 

{0; 0; 0} 


PQ 

{—i; —i ± j) 

{-1} 


(I-4QR) 1 / 2 

{— l; — i ±j) 

{1; 1 ± V3} 


I-QP 

{— i; —i ± j) 

{0; 0; 0} 


QP 

{— i; —i ± j) 

{-1} 


Fact 5, along with the preceding four facts, permits a complete classification of the time-invariant 
linear multivariable unity feedback system examined in this paper. These facts imply that 1) if the given 
plant P and the compensator C are unstable and they have common unstable poles, then to guarantee the 
stability of the feedback system it is required to ascertain that (i) both the closed-loop transfer matrix from 
the command input ul to the compensator output yl (which is Q) and the closed-loop transfer matrix 
from the disturbance input u2 to the plant output y2 (which is R) are stable, and (ii) both 
(I - 4RQ) 1 / 2 and (I - 4QR) 1 / 2 are stable (thus guaranteeing that the closed- loop transfer matrices from 
the command input to the plant output and from the disturbance input to the compensator output are 
stable); 2) if the plant and the compensator are unstable but they do not have common unstable poles, then 
it is required to ascertain only the stability of the associated Q and R; 3) if the plant is unstable but the 
compensator is stable, then it is required to only ascertain the stability of the associated R (the associated 
Q will be stable automatically); 4) if the plant is stable but the compensator is unstable, then we only need 
to make sure that the associated Q is stable (the stability of the associated R will be automatically 
satisfied); 5) if the given plant and the compensator are both stable, then we only need to ascertain that 
either the associated R or Q is stable (both will be stable automatically). These facts are summarized in 
table 4. For the sake of completeness, the SISO case is also included in the table. 
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TABLE 4.- CLASSIFICATION OF MIMO UNITY FEEDBACK SYSTEMS USING Q- AND 

R-PARAMETRIZATION. 


Stability characteristics 
of the plant P and the 
compensator C 

Necessary and sufficient conditions for the 
feedback system to be stable 


SISO 

MIMO 

P: unstable 
C: unstable 

CandP have common 
unstable poles 

Q = C(I + PC) -1 stable 
and 

R = P(I + CP) -1 stable 

Q and R stable, and 
(I-4RQ)l/2 and 
(I-4QR)l/2 stable 

P: unstable 
C: unstable 

CandP do not have 
common unstable poles 

QandR stable 

QandR stable 

P: unstable 
C: stable 

R stable 

R stable 

P: stable 
C: unstable 

Q stable 

Q stable 

P: stable 
C: stable 

Q stable or 
R stable 

Q stable or 
R stable 


4. COMPARISON WITH STABILITY TESTS USING THE CHARACTERISTIC 

POLYNOMIAL 


Using a minimal state-space realization, it has been known for many years (Hsu et al. (1968)) that 
the characteristic polynomial of the feedback system, A(s) is given by 

A(s) = Ac(s) Ap(s)det[I + PC] (24) 

where Ac(s) and Ap(s) are respectively the characteristic polynomial of the compensator and the plant. 
Note that det[I + PC] = det[I + CP]. Table 5 shows the results of applying equation (24) to the five 
examples listed in table 1. As expected, the results concerning the stability of the feedback system are 
identical to those shown in table 2 discussed earlier. 
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TABLE 5.- CHARACTERISTIC POLYNOMIAL OF THE 
FIVE EXAMPLES IN TABLE 1. 


Examples in 
table 1 

Characteristic polynomial of the 
feedback system 

#1 

(s + 1)2( S 2 + s + 2) 

#2 

s(s 2 + 3s + l)(s 2 + s + 1) 

#3 

s(s 2 + 3s + l)(s 2 + s + 1) 

#4 

(s - l)2(s 2 + 2s + 5)2 

#5 

(s + 1) 2 ( S 2 + 2s + 2) 


Notice that, to apply equation (24), a minimal state-space realization (or equivalently a co-prime 
factorization (see Desoer et al. (1975)) of the plant and the compensator are required. By contrast, how- 
ever, no such procedure is required in table 4 with Q- and R-Parametrization. 


5. ILLUSTRATIVE EXAMPLES 


The following simple examples serve to illustrate some of the applications that make use of the 
ideas shown in table 4. 

Example 1 . To begin with, consider a very simple unstable SISO plant as shown in figure 4. 

COMPENSATOR PLANT 


ul el 

^ 

C(s) 

yl 



u2 

k ^ 

s 



e2 

s-1 


Figure 4 - A simple SISO feedback system. 

Can the compensators of the forms, (a) C(s) = Cl(s) = k(s + l)/s and (b) C(s) = C2(s) = 
k(s + l)/(s -0.5), where k being a constant, stabilize the plant? Generate all the compensators that 
stabilize the plant 

We note that since both the plant and the compensators (a) and (b) are unstable for this SISO feed- 
back system, we need, according to table 4, to ascertain that both the associated Q and R are stable. For 
Cl(s), we found that die associated Q1 and Rl, equations (3) and (5), cannot be both stable for any value 
of k. We conclude therefore that Cl(s) cannot stabilize the plant. For C2(s), the associated Q2 and R2 
are both stable with values of k > 1.5. Thus, the plant can be stabilized by the C2(s) with k > 1.5. 
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To generate all the compensators, stable or otherwise by themselves, which stabilize the unstable 
SISO plant, we must again ascertain that both the associated Q and R are stable. To meet this require- 
ment, we first choose a general, stable, and proper Q1 such as 

q. _ pi(s) _ pOs* + pls i_1 + p2s*~ 2 4 + pi 

qi(s) s* + qls i_1 + q2s* -2 + • • • +qi 

with qi(s) being Hurwitz. We then enforce the stability requirement of the associated Ri, which from 
equation (12) is 


R . _ n(s)[d(s)qi(s) - n(s)pi(s)] 

[d(s)] 2 qi(s) 

where the plant P(s) is denoted by P(s) = n(s)/d(s). This permits us to express the coefficients of pi(s) 
in terms of the coefficients of qi(s). Finally the associated compensators, Ci, that stabilize the plant are 
obtained. 


_ d(s)pi(s) 
d(s)qi(s) - n(s)pi(s) 

The following compensators are generated from the above process with i = 1 and i = 2. 
a) i = 1: 


Q1= £ s + q!p ; ql > 0, pi = -(1 + ql), pO = 1 + ql 


~<\ l 


+ ql 


Cl =-(l + 1/ql) 


b) i = 2: 

Q2 = P 0s +pls + p2 . qi > o, q2 > 0, pO arbitrary, 
s 2 + qls + q2 

pi = (1 + ql + q2) - 2p0, and 

p2 = —(1 + ql + q2) + pO 

R2 (1 - P 0)s[s-q2/(l -pQ>] 
s 2 + qls + q2 

C2 _ p0s 2 + pls + p2 pO s - [1 - (1 + ql + q2)/p0] 

(l-p0)(s-l)[s-q2/(l-p0)] 1-pO s-q2/(l - pO) 

The above 3-parameter compensators can be reduced to various 2-parameter first-order compensators with 
a proper selection of pO, pi, and p2. 


i 
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1) pO = 0: 


C2 = 


1 + ql + q2 
s -q2 


2) p2 = 0: 


3) pi = 0: 


Q 2 _ 1 + ql + q2 


ql + q2 s + q2/(ql + q2) 


po _ 1 + ql + q2 s + 1 

~ 1 -ql -q2 s + 2q2/(ql + q2 - 1) 

In particular, if it is now restricted to 

ql = 1 - 5q2 , q2 < 0.2 

then the family of compensators C2(s) of (b) is produced. 

02 = k s + - k > 1 5 
C s-0.5’ D 


Example 2 . Consider a simple unstable MIMO plant P(s) and the compensators of the proportional plus 
integral form: 


■ 1 1 ■ 


, , k2 



kl + — 0 

s + 1 s- 1 


s 

1 

c = 

k4 

0 — l — 


0 k3 + — 

S - 1 


s 


Can C stabilize P? 

Both the C and P are unstable, but they do not have common unstable poles. Thus from table 4, 
C stabilizes P iff QandR are stable. We found that both QandR are stable if 

kl > -1 , k2 > 0 and k3 > 1 , k4 > 0 

Therefore C can stabilize P if ki (i = 1 , 4) are satisfied by the above conditions. In particular, if we set 
kl = k2 = 1 and k3 = k4 = 2, then we reduce the compensators to the special one shown in the 
example #1 in table 1 discussed previously. 

Example 3 . Given a two-input, two-output plant P(s), 


P(s) = 


s - 2 

(s-lXs + 2) 
3(s - 2) 

(s - l)(s + 2) 


2(s - 2) 

(s - l)(s + 2) 
4(s-2) 
(S'- l)(s + 2) 


Can the plant be stabilized with simple proportional compensators of the following form? 
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C(s) = 


kl 0‘ 
0 k2 


Since the plant is unstable but the compensators are stable, we need, according to table 4, only to ascertain 
the stability of the associated R. We found that R is 

s - 2 r (s + 2)(s - 1) - 2(s - 2)k2 2(s + 2)(s-l) 

r(s) [ 3(s + 2)(s - 1) 4(s + 2)(s — 1) - 2(s - 2)kl 

r(s) = s 4 + (a + b)s 3 + (ab - 2a - 2b - 6klk2)s 2 + 4(6klk2 - ab)s - 4(6klk2 - ab) 


where a = 1 + kl, and b = 1 + 4k2. We see from the last two terms of r(s) that r(s) is not Hurwitz 
and r(s) cannot be the product of a Hurwitz and (s - 2). Thus R is always unstable and hence C(s) 
cannot stabilize the plant. 


Example 4 . Generate all the plants which are stabilizable by the following family of compensators, 


C(s) = 


nls + nO 
s + dO 

0 


0 

n3s + n2 » 
s + dl 


dO > 0, dl > 0 


We recall that this is the sort of problem which provides a motivation behind our introduction of the 
R-Parametrization described earlier in the paper. The family of compensators is stable by itself. There- 
fore, the feedback system is stable iff R is stable. Let us therefore choose a sequence of stable and 
proper Rs such as 

Rs = {R0, Rl, R2, . . ., Ri, . . .} 

where 

Ri = -r|— (kOs* + Kls i_1 + K2s i_2 + • • ■ + Ki) 
n(s) ' 1 

and where K0, Kl, K2, . . ., Ki are arbitrary 2x2 constant matrices, and ri(s) is a Hurwitz polynomial 
of order i in s, i.e.. 


ri(s) = s 1 + rls 1 1 + r2s‘ 2 + • • • + ri 

The sequence of the plants, Ps, associated with the compensator C and the Rs is obtained. 

Ps = {P0, P1,P2, ..., Pi,...} 


where 


Pi = Ri(I - CRi) -1 

Since Rs is a general and proper sequence, Ps is a sequence consisting of all the plants stabilizable by the 
given stable compensator. The following plants are generated from the above process with i = 0, and 
i = 1. 
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1) i = 0: 


R0 = K0 = 


kl k2 
k3k4j’ 


ki are arbitrary constants 


PO = RO(I - CRO)" 1 

1 [ tOl t02 

aO(s) t03 t04 

where 

aO(s) = [(1 - k4n3)s + (dl -k4n2)][(l -klnl)s + (dO-klnO)] -k2k3(nls + nO)(n3s + n2) 
tOl = (s + dO){[kl(l - k4n3) + k2k3n3]s + kl(dl -k4n2) + k2k3n2} 
t02 = (s + dl){ [k2(l - klnl) + klk2nl]s + k2(d0 - klnO) + klk2n0} 
t03 = (s + dO){ [k3(l - k4n3) + k3k4n3]s + k3(dl - k4n2) + k3k4n2} 
and 

t04 = (s + dl){[k4(l - klnl) + k2k3nl]s + k4(d0 - klnO) + k2k3n0) 

2) i = 1: 

R1 = + ^ (KOs + Kl) , Kl = ^ , ki are arbitrary constants 

PI =R1(I-CR1)" 1 

1 [ tH tl2 ' 

al(s) [ tl3 tl4 . 

where 

al(s) = [(s + rl)(s + dl) - (n3s + n2)(k4s + k8)][(s + rl)(s + dO) - (nls + nO)(kls + k5)] 

- (k2s + k6)(k3s + k7)(nls + n0)(n3s + n2) 

til = (s + dO){(kls + k5)[(s + rl)(s + dl) - (n3s + n2)(k4s + k8)] 

+ (k2s + k6)(n3s + n2)(k3s + k7)} 

tl2 = (s + dl){(kls + k5)(nls + n0)(k2s + k6) + (k2s + k6)[(s + rl)(s + dO) 

- (nls + nO)(kls + k5)]} 

tl3 = (s + d0){(k3s + k7)[(s + rl)(s + dl) - (n3s + n2)(k4s + k8)] 

+ (k4s + k8)(n3s + n2)(k3s + k7)} 
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and 


tl4 = (s + dl){(k3s + k7)(nls + n0)(k2s + k6) + (k4s + k8)[(s + rl)(s + dO) 

- (nls + nO)(kls + k5)]} 

We note in passing that, for nO = dO = n2 = dl = 0, the compensators reduce to the special proportional 
type of compensators of example 3. As expected, under this condition, the plant in example 3 is not a 
member of Ps. However, the following plant, 


P(s) = 


(s-4) 

(s- l)(s + 4) L 


1 2 
3 4 


for instance, is stabilizable by the special proportional type of compensators, and we find that this plant is 
indeed a member of Ps. 


6. CONCLUSIONS 


A classification of all the compensators which stabilize a given plant, stable or unstable, in a linear, 
time-invariant, MIMO unity feedback system is developed. This classification is made possible by intro- 
ducing a new parametrization referred herein as R-Parametrization, which is dual of the familiar 
Q-Parametrization. It is shown that (i) if the plant and the compensators are both unstable and they 
have common unstable poles, then the feedback system is stable iff (a) Q and R are stable, and 
(b) (I - 4RQ)!/ 2 and (I - 4QR) 1 / 2 are stable, (ii) if the plant and the compensators are both unstable 
but they do not have common unstable poles, then the feedback system is stable iff both Q and R are 
stable, (iii) if the plant is unstable but the compensators are stable, then the feedback system is stable iff R 
is stable, (iv) if the plant is stable but the compensators are unstable, then the feedback system is stable iff 
Q is stable, and (v) if the plant and the compensators are both stable, then the feedback system is stable 
iff either Q or R is stable. 
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